Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets.

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

Genome-wide study investigating effector genes and polygenic prediction for kidney function in persons with ancestry from Africa and the Americas.

Chronic kidney disease is a leading cause of death and disability globally and impacts individuals of African ancestry (AFR) or with ancestry in the Americas (AMS) who are under-represented in genome-wide association studies (GWASs) of kidney function. To address this bias, we conducted a large meta-analysis of GWASs of estimated glomerular filtration rate (eGFR) in 145,732 AFR and AMS individuals. We identified 41 loci at genome-wide significance (p < 5 × 10-8), of which two have not been previously reported in any ancestry group. We integrated fine-mapped loci with epigenomic and transcriptomic resources to highlight potential effector genes relevant to kidney physiology and disease, and reveal key regulatory elements and pathways involved in renal function and development. We demonstrate the varying but increased predictive power offered by a multi-ancestry polygenic score for eGFR and highlight the importance of population diversity in GWASs and multi-omics resources to enhance opportunities for clinical translation for all.

Towards Faster Gene Expression Prediction via Dimensionality Reduction and Feature Selection.

The majority of genes have a genetic component to their expression. Elastic nets have been shown effective at predicting tissue-specific, individual-level gene expression from genotype data. We apply principal component analysis (PCA), linkage disequilibrium pruning, or the combination of the two to reduce, or generate, a lower-dimensional representation of the genetic variants used as inputs to the elastic net models for the prediction of gene expression. Our results show that, in general, elastic nets attain their best performance when all genetic variants are included as inputs; however, a relatively low number of principal components can effectively summarize the majority of genetic variation while reducing the overall computation time. Specifically, 100 principal components reduce the computational time of the models by over 80% with only an 8% loss in R2. Finally, linkage disequilibrium pruning does not effectively reduce the genetic variants for predicting gene expression. As predictive models are commonly made for over 27,000 genes for more than 50 tissues, PCA may provide an effective method for reducing the computational burden of gene expression analysis.

Effects of iron substitution and anti-site disorder on crystal structures, vibrational, optical and magnetic properties of double perovskites Sr2(Fe1-xNix)TeO6.

The double-perovskite series, Sr2(Fe1-xNix)TeO6 (x = 0, 0.25, 0.50, 0.75, and 1) has been synthesized in polycrystalline form by solid-state reaction at 1300 K in air. Their crystal structures were probed by powder X-ray diffraction at room temperature. Rietveld analysis revealed that all samples crystallize in the monoclinic space group I2/m. The double-perovskite structures ideally contain two alternating types of octahedra (Fe/Ni)2dO6 and (Te)2aO6, tilted in the system (a-a-c0). However, the refinements have shown a complex distribution of all three cations over the two available octahedral sites; 2d (½, ½, 0) and 2a (0, 0, 0). Raman spectroscopy further complements the obtained results, by revealing a tiny increase of the wavenumber of some Raman modes when Fe is substituted by Ni. The optical characteristics of the series were determined by fitting diffuse reflectance UV/Vis spectra enabling the optical band gaps to be derived from Tauc method and derivation of absorption spectra fitting (DASF) techniques. Analyses of the obtained 57Fe Mössbauer hyperfine parameters at room temperature of samples with compositions x = 0, 0.25, 0.50 and 0.75 reveal the presence of Fe3+ in high-spin state with an anti-site disorder of Fe-Ni-Te cations in distorted octahedral environments (site 2d and 2a). The results show that significant correlations exist between the crystal structures and physical properties of double perovskites containing B site transition elements of different charge and size. Temperature-dependent magnetic susceptibility data show magnetic transitions below 40(1) K (38(1) K, 31(1) K, 25(1) K, 20(1) K, and 35(1) K for x = 0, 0.25, 0.50, 0.75, and 1, respectively. A divergence between FC and ZFC curves for all compositions has been observed. The results show that the ground states of the doped materials might be spin glasses or magnetically ordered.

Investigation of Copy Number Variation in South African Patients With Congenital Heart Defects.

BACKGROUND: Congenital heart disease (CHD) is a leading non-infectious cause of pediatric morbidity and mortality worldwide. Although the etiology of CHD is poorly understood, genetic factors including copy number variants (CNVs) contribute to the risk of CHD in individuals of European ancestry. The presence of rare CNVs in African CHD populations is unknown. This study aimed to identify pathogenic and likely pathogenic CNVs in South African patients with CHD. METHODS: Genotyping was performed on 90 patients with nonsyndromic CHD using the Affymetrix CytoScan HD platform. These data were used to identify large, rare CNVs in known CHD-associated genes and candidate genes. RESULTS: We identified eight CNVs overlapping known CHD-associated genes (GATA4, CRKL, TBX1, FLT4, B3GAT3, NSD1) in six patients. The analysis also revealed CNVs encompassing five candidate genes likely to play a role in the development of CHD (DGCR8, KDM2A, JARID2, FSTL1, CYFIP1) in five patients. One patient was found to have 47, XXY karyotype. We report a total discovery yield of 6.7%, with 5.6% of the cohort carrying pathogenic or likely pathogenic CNVs expected to cause the observed phenotypes. CONCLUSIONS: In this study, we show that chromosomal microarray is an effective technique for identifying CNVs in African patients diagnosed with CHD and have demonstrated results similar to previous CHD genetic studies in Europeans. Novel potential CHD genes were also identified, indicating the value of genetic studies of CHD in ancestrally diverse populations.

Telomere therapy for chronic kidney disease.

Chronic kidney disease (CKD) is estimated to affect almost 10% of individuals worldwide and is one of the leading causes of morbidity and mortality. Renal fibrosis, a central pathway in CKD progression (irrespective of etiology), is associated with shortened or dysfunctional telomeres in animal studies. Telomeres are specialized nucleoprotein structures located at the chromosome end that maintain genomic integrity. The mechanisms of associations between telomere length and CKD have not yet been fully elucidated, however, CKD patients with shorter telomere length may have decreased renal function and a higher mortality rate. A plethora of ongoing research has focused on possible therapeutic applications of telomeres with the overall goal to preserve telomere length as a therapy to treat CKD.

Chronic kidney disease or CKD is one of the leading causes of illness and death worldwide. Scarring of the kidney tissue that occurs in CKD has been associated with shorter telomeres in studies using rats. Telomeres, said to act as the cellular 'shoelace caps', maintain the structure of chromosomes, allowing for genetic material inside cells to divide correctly. The length of telomeres (TL) is influenced by diverse factors such as genetics and lifestyle. The underlying processes for the associations between TL and CKD are still not understood, however, patients with CKD and shorter TL have reduced kidney function and an increased death rate. Therefore, research is focused on possible ways to preserve TL and treat CKD.

Adapting Random Forests to Predict Obesity-Associated Gene Expression.

Random forests (RFs) are effective at predicting gene expression from genotype data. However, a comparison of RF regressors and classifiers, including feature selection and encoding, has been under-explored in the context of gene expression prediction. Specifically, we examine the role of ordinal or one-hot encoding and of data balancing via oversam-pling in the prediction of obesity-associated gene expression. Our work shows that RFs compete with PrediXcan in the prediction of obesity-associated gene expression in subcutaneous adipose tissue, a highly relevant tissue to obesity. Additionally, RFs generate predictions for obesity-associated genes where PrediXcan fails to do so.

Kidney omics in hypertension: from statistical associations to biological mechanisms and clinical applications.

Hypertension is a major cardiovascular disease risk factor and contributor to premature death globally. Family-based investigations confirmed a significant heritable component of blood pressure (BP), whereas genome-wide association studies revealed >1000 common and rare genetic variants associated with BP and/or hypertension. The kidney is not only an organ of key relevance to BP regulation and the development of hypertension, but it also acts as the tissue mediator of genetic predisposition to hypertension. The identity of kidney genes, pathways, and related mechanisms underlying the genetic associations with BP has started to emerge through integration of genomics with kidney transcriptomics, epigenomics, and other omics as well as through applications of causal inference, such as Mendelian randomization. Single-cell methods further enabled mapping of BP-associated kidney genes to cell types, and in conjunction with other omics, started to illuminate the biological mechanisms underpinning associations of BP-associated genetic variants and kidney genes. Polygenic risk scores derived from genome-wide association studies and refined on kidney omics hold the promise of enhanced diagnostic prediction, whereas kidney omics-informed drug discovery is likely to contribute new therapeutic opportunities for hypertension and hypertension-mediated kidney damage.

Contributions of obesity to kidney health and disease: insights from Mendelian randomization and the human kidney transcriptomics.

AIMS: Obesity and kidney diseases are common complex disorders with an increasing clinical and economic impact on healthcare around the globe. Our objective was to examine if modifiable anthropometric obesity indices show putatively causal association with kidney health and disease and highlight biological mechanisms of potential relevance to the association between obesity and the kidney. METHODS AND RESULTS: We performed observational, one-sample, two-sample Mendelian randomization (MR) and multivariable MR studies in ∼300 000 participants of white-British ancestry from UK Biobank and participants of predominantly European ancestry from genome-wide association studies. The MR analyses revealed that increasing values of genetically predicted body mass index and waist circumference were causally associated with biochemical indices of renal function, kidney health index (a composite renal outcome derived from blood biochemistry, urine analysis, and International Classification of Disease-based kidney disease diagnoses), and both acute and chronic kidney diseases of different aetiologies including hypertensive renal disease and diabetic nephropathy. Approximately 13-16% and 21-26% of the potentially causal effect of obesity indices on kidney health were mediated by blood pressure and type 2 diabetes, respectively. A total of 61 pathways mapping primarily onto transcriptional/translational regulation, innate and adaptive immunity, and extracellular matrix and metabolism were associated with obesity measures in gene set enrichment analysis in up to 467 kidney transcriptomes. CONCLUSIONS: Our data show that a putatively causal association of obesity with renal health is largely independent of blood pressure and type 2 diabetes and uncover the signatures of obesity on the transcriptome of human kidney.

Polygenic basis and biomedical consequences of telomere length variation.

Telomeres, the end fragments of chromosomes, play key roles in cellular proliferation and senescence. Here we characterize the genetic architecture of naturally occurring variation in leukocyte telomere length (LTL) and identify causal links between LTL and biomedical phenotypes in 472,174 well-characterized UK Biobank participants. We identified 197 independent sentinel variants associated with LTL at 138 genomic loci (108 new). Genetically determined differences in LTL were associated with multiple biological traits, ranging from height to bone marrow function, as well as several diseases spanning neoplastic, vascular and inflammatory pathologies. Finally, we estimated that, at the age of 40 years, people with an LTL >1 s.d. shorter than the population mean had a 2.5-year-lower life expectancy compared with the group with ≥1 s.d. longer LDL. Overall, we furnish new insights into the genetic regulation of LTL, reveal wide-ranging influences of LTL on physiological traits, diseases and longevity, and provide a powerful resource available to the global research community.

Plasma Proteomics of Renal Function: A Transethnic Meta-Analysis and Mendelian Randomization Study.

BACKGROUND: Studies on the relationship between renal function and the human plasma proteome have identified several potential biomarkers. However, investigations have been conducted largely in European populations, and causality of the associations between plasma proteins and kidney function has never been addressed. METHODS: A cross-sectional study of 993 plasma proteins among 2882 participants in four studies of European and admixed ancestries (KORA, INTERVAL, HUNT, QMDiab) identified transethnic associations between eGFR/CKD and proteomic biomarkers. For the replicated associations, two-sample bidirectional Mendelian randomization (MR) was used to investigate potential causal relationships. Publicly available datasets and transcriptomic data from independent studies were used to examine the association between gene expression in kidney tissue and eGFR. RESULTS: In total, 57 plasma proteins were associated with eGFR, including one novel protein. Of these, 23 were additionally associated with CKD. The strongest inferred causal effect was the positive effect of eGFR on testican-2, in line with the known biological role of this protein and the expression of its protein-coding gene (SPOCK2) in renal tissue. We also observed suggestive evidence of an effect of melanoma inhibitory activity (MIA), carbonic anhydrase III, and cystatin-M on eGFR. CONCLUSIONS: In a discovery-replication setting, we identified 57 proteins transethnically associated with eGFR. The revealed causal relationships are an important stepping stone in establishing testican-2 as a clinically relevant physiological marker of kidney disease progression, and point to additional proteins warranting further investigation.

Epigenome-wide association study of kidney function identifies trans-ethnic and ethnic-specific loci.

BACKGROUND: DNA methylation (DNAm) is associated with gene regulation and estimated glomerular filtration rate (eGFR), a measure of kidney function. Decreased eGFR is more common among US Hispanics and African Americans. The causes for this are poorly understood. We aimed to identify trans-ethnic and ethnic-specific differentially methylated positions (DMPs) associated with eGFR using an agnostic, genome-wide approach. METHODS: The study included up to 5428 participants from multi-ethnic studies for discovery and 8109 participants for replication. We tested the associations between whole blood DNAm and eGFR using beta values from Illumina 450K or EPIC arrays. Ethnicity-stratified analyses were performed using linear mixed models adjusting for age, sex, smoking, and study-specific and technical variables. Summary results were meta-analyzed within and across ethnicities. Findings were assessed using integrative epigenomics methods and pathway analyses. RESULTS: We identified 93 DMPs associated with eGFR at an FDR of 0.05 and replicated 13 and 1 DMPs across independent samples in trans-ethnic and African American meta-analyses, respectively. The study also validated 6 previously published DMPs. Identified DMPs showed significant overlap enrichment with DNase I hypersensitive sites in kidney tissue, sites associated with the expression of proximal genes, and transcription factor motifs and pathways associated with kidney tissue and kidney development. CONCLUSIONS: We uncovered trans-ethnic and ethnic-specific DMPs associated with eGFR, including DMPs enriched in regulatory elements in kidney tissue and pathways related to kidney development. These findings shed light on epigenetic mechanisms associated with kidney function, bridging the gap between population-specific eGFR-associated DNAm and tissue-specific regulatory context.

Uncovering genetic mechanisms of hypertension through multi-omic analysis of the kidney.

The kidney is an organ of key relevance to blood pressure (BP) regulation, hypertension and antihypertensive treatment. However, genetically mediated renal mechanisms underlying susceptibility to hypertension remain poorly understood. We integrated genotype, gene expression, alternative splicing and DNA methylation profiles of up to 430 human kidneys to characterize the effects of BP index variants from genome-wide association studies (GWASs) on renal transcriptome and epigenome. We uncovered kidney targets for 479 (58.3%) BP-GWAS variants and paired 49 BP-GWAS kidney genes with 210 licensed drugs. Our colocalization and Mendelian randomization analyses identified 179 unique kidney genes with evidence of putatively causal effects on BP. Through Mendelian randomization, we also uncovered effects of BP on renal outcomes commonly affecting patients with hypertension. Collectively, our studies identified genetic variants, kidney genes, molecular mechanisms and biological pathways of key relevance to the genetic regulation of BP and inherited susceptibility to hypertension.

Nanoparticle-Enabled Enrichment of Longitudinal Blood Proteomic Fingerprints in Alzheimer's Disease.

Blood-circulating biomarkers have the potential to detect Alzheimer's disease (AD) pathology before clinical symptoms emerge and to improve the outcomes of clinical trials for disease-modifying therapies. Despite recent advances in understanding concomitant systemic abnormalities, there are currently no validated or clinically used blood-based biomarkers for AD. The extremely low concentration of neurodegeneration-associated proteins in blood necessitates the development of analytical platforms to address the "signal-to-noise" issue and to allow an in-depth analysis of the plasma proteome. Here, we aimed to discover and longitudinally track alterations of the blood proteome in a transgenic mouse model of AD, using a nanoparticle-based proteomics enrichment approach. We employed blood-circulating, lipid-based nanoparticles to extract, analyze and monitor AD-specific protein signatures and to systemically uncover molecular pathways associated with AD progression. Our data revealed the existence of multiple proteomic signals in blood, indicative of the asymptomatic stages of AD. Comprehensive analysis of the nanoparticle-recovered blood proteome by label-free liquid chromatography-tandem mass spectrometry resulted in the discovery of AD-monitoring signatures that could discriminate the asymptomatic phase from amyloidopathy and cognitive deterioration. While the majority of differentially abundant plasma proteins were found to be upregulated at the initial asymptomatic stages, the abundance of these molecules was significantly reduced as a result of amyloidosis, suggesting a disease-stage-dependent fluctuation of the AD-specific blood proteome. The potential use of the proposed nano-omics approach to uncover information in the blood that is directly associated with brain neurodegeneration was further exemplified by the recovery of focal adhesion cascade proteins. We herein propose the integration of nanotechnology with already existing proteomic analytical tools in order to enrich the identification of blood-circulating signals of neurodegeneration, reinvigorating the potential clinical utility of the blood proteome at predicting the onset and kinetics of the AD progression trajectory.

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals.

Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to ~1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency ≤ 0.01) variant BP associations (P < 5 × 10-8), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were ~8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets.

Hypertension and renin-angiotensin system blockers are not associated with expression of angiotensin-converting enzyme 2 (ACE2) in the kidney.

AIMS: Angiotensin-converting enzyme 2 (ACE2) is the cellular entry point for severe acute respiratory syndrome coronavirus (SARS-CoV-2)-the cause of coronavirus disease 2019 (COVID-19). However, the effect of renin-angiotensin system (RAS)-inhibition on ACE2 expression in human tissues of key relevance to blood pressure regulation and COVID-19 infection has not previously been reported. METHODS AND RESULTS: We examined how hypertension, its major metabolic co-phenotypes, and antihypertensive medications relate to ACE2 renal expression using information from up to 436 patients whose kidney transcriptomes were characterized by RNA-sequencing. We further validated some of the key observations in other human tissues and/or a controlled experimental model. Our data reveal increasing expression of ACE2 with age in both human lungs and the kidney. We show no association between renal expression of ACE2 and either hypertension or common types of RAS inhibiting drugs. We demonstrate that renal abundance of ACE2 is positively associated with a biochemical index of kidney function and show a strong enrichment for genes responsible for kidney health and disease in ACE2 co-expression analysis. CONCLUSION: Our results indicate that neither hypertension nor antihypertensive treatment is likely to alter the expression of the key entry receptor for SARS-CoV-2 in the human kidney. Our data further suggest that in the absence of SARS-CoV-2 infection, kidney ACE2 is most likely nephro-protective but the age-related increase in its expression within lungs and kidneys may be relevant to the risk of SARS-CoV-2 infection.

Human Y Chromosome Exerts Pleiotropic Effects on Susceptibility to Atherosclerosis.

OBJECTIVE: The male-specific region of the Y chromosome (MSY) remains one of the most unexplored regions of the genome. We sought to examine how the genetic variants of the MSY influence male susceptibility to coronary artery disease (CAD) and atherosclerosis. Approach and Results: Analysis of 129 133 men from UK Biobank revealed that only one of 7 common MSY haplogroups (haplogroup I1) was associated with CAD-carriers of haplogroup I1 had ≈11% increase in risk of CAD when compared with all other haplogroups combined (odds ratio, 1.11; 95% CI, 1.04-1.18; P=6.8×10-4). Targeted MSY sequencing uncovered 235 variants exclusive to this haplogroup. The haplogroup I1-specific variants showed 2.45- and 1.56-fold respective enrichment for promoter and enhancer chromatin states, in cells/tissues relevant to atherosclerosis, when compared with other MSY variants. Gene set enrichment analysis in CAD-relevant tissues showed that haplogroup I1 was associated with changes in pathways responsible for early and late stages of atherosclerosis development including defence against pathogens, immunity, oxidative phosphorylation, mitochondrial respiration, lipids, coagulation, and extracellular matrix remodeling. UTY was the only Y chromosome gene whose blood expression was associated with haplogroup I1. Experimental reduction of UTY expression in macrophages led to changes in expression of 59 pathways (28 of which overlapped with those associated with haplogroup I1) and a significant reduction in the immune costimulatory signal. CONCLUSIONS: Haplogroup I1 is enriched for regulatory chromatin variants in numerous cells of relevance to CAD and increases cardiovascular risk through proatherosclerotic reprogramming of the transcriptome, partly through UTY.

Uncovering genetic mechanisms of kidney aging through transcriptomics, genomics, and epigenomics.

Nephrons scar and involute during aging, increasing the risk of chronic kidney disease. Little is known, however, about genetic mechanisms of kidney aging. We sought to define the signatures of age on the renal transcriptome using 563 human kidneys. The initial discovery analysis of 260 kidney transcriptomes from the TRANScriptome of renaL humAn TissuE Study (TRANSLATE) and the Cancer Genome Atlas identified 37 age-associated genes. For 19 of those genes, the association with age was replicated in 303 kidney transcriptomes from the Nephroseq resource. Surveying 42 nonrenal tissues from the Genotype-Tissue Expression project revealed that, for approximately a fifth of the replicated genes, the association with age was kidney-specific. Seventy-three percent of the replicated genes were associated with functional or histological parameters of age-related decline in kidney health, including glomerular filtration rate, glomerulosclerosis, interstitial fibrosis, tubular atrophy, and arterial narrowing. Common genetic variants in four of the age-related genes, namely LYG1, PPP1R3C, LTF and TSPYL5, correlated with the trajectory of age-related changes in their renal expression. Integrative analysis of genomic, epigenomic, and transcriptomic information revealed that the observed age-related decline in renal TSPYL5 expression was determined both genetically and epigenetically. Thus, this study revealed robust molecular signatures of the aging kidney and new regulatory mechanisms of age-related change in the kidney transcriptome.

Trans-ethnic kidney function association study reveals putative causal genes and effects on kidney-specific disease aetiologies.

Chronic kidney disease (CKD) affects ~10% of the global population, with considerable ethnic differences in prevalence and aetiology. We assemble genome-wide association studies of estimated glomerular filtration rate (eGFR), a measure of kidney function that defines CKD, in 312,468 individuals of diverse ancestry. We identify 127 distinct association signals with homogeneous effects on eGFR across ancestries and enrichment in genomic annotations including kidney-specific histone modifications. Fine-mapping reveals 40 high-confidence variants driving eGFR associations and highlights putative causal genes with cell-type specific expression in glomerulus, and in proximal and distal nephron. Mendelian randomisation supports causal effects of eGFR on overall and cause-specific CKD, kidney stone formation, diastolic blood pressure and hypertension. These results define novel molecular mechanisms and putative causal genes for eGFR, offering insight into clinical outcomes and routes to CKD treatment development.